Fragility of String Orders 
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One-dimensional gapped systems are often characterized by a 'hidden' non-local order parameter, 
the so-called string order. Due to the gap, thermodynamic properties are robust against a weak 
higher-dimensional coupling between such chains or ladders. To the contrary, we find that the string 
order is not stable and decays for arbitrary weak inter-chain or inter-ladder coupling. We investigate 
the vanishing of the order for three different systems: spin-one Haldane chains, band insulators, and 
the transverse field Ising model. Using perturbation theory and bosonization, we show that the 
fragility of the string order arises from non-local commutation relations between the non-local order 
parameter and the perturbation. 
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Since the pioneering work of Landau— that led to the 
theory of 2nd order phase transitions, the concept of or- 
der parameter has become one of the main paradigms 
in condensed matter physics. The central observation is 
that when a symmetry is spontaneously broken and long 
range order appears into the system, certain correlation 
functions of operators^. A{x) do not decay with distance. 



order" defined by the non-local correlation function 



lim (A(x)A(y)) = c ^ 0. 

|x-y|^oo 



(1) 



and it is possible to define an order parameter as (A) — 
^Jc. For example, it is known that a classical three- 
dimensional ferromagnet undergoes a phase transition, 
becoming ordered (all the spins point in the same direc- 
tion) below a certain critical temperature. In this case, 
the order parameter can be identified with the sponta- 
neous magnetization (S(a;)) — m and the spin-spin corre- 
lation function {S{x)S{y)) — is finite also in the limit 
\x — y\ ^ oo. This simple way to describe and discrimi- 
nate the different phases of matter has been successfully 
applied to a plethora of different systems both at finite 
and zero temperature. Classical and quantum magnetism 
(with many possible order parameters: uniform and stag- 
gered magnetization, dimerization...), superconductivity 
(Cooper pair amplitude), superfluidity (condensate am- 
plitude) are only few of the many successful applications 
of the Landau paradigm. 

In the context of low dimensional quantum systems, it 
has been also useful, under some circumstances, to iden- 
tify 'non-local' order parameters which cannot be written 
in the form ([T]) . Typically, this is the case of some gapped 
one-dimensional Hamiltonians with no local symmetries 
that are spontaneously broken and with any two-point 
correlation function that decays exponentially. In such 
systems an hidden long range order can nevertheless be 
present and this is encoded in the long distance behav- 
ior of certain non-local operators. Notice that a 'non- 
local' order parameter is not directly accessible to any 
experimental probe but can be equally used to mark the- 
oretically the boundaries of the different phases. In this 
context, an important example is the so-called "string- 
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|x-a|-n 



(A(x) n i?(z) A{y)) ^ 0, 
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where the operator Hzes ^(^) ^-cts on a line (the string 
•Sx.y) connecting the points x and y. One dimensional 
examples of this family (discussed further below) includes 
the spin-1 chain^ and spin- 1/2 ladders^, the transverse- 
field quantum Ising chain (dual order), band- and Mott- 
Fermionic insulators^ and, among Bosonic systems, some 
parameter regime of the Bose-Hubbard Hamiltonian^. 

Another important sub-class of these more exotic types 
of systems is characterized by the so-called 'topological 
order' (e.g. fractional quantum Hall fiuids). This concept 
can be defined by the ground state degeneracy on non- 
trivial manifolda^i^i^. Here we will only consider systems, 
as listed above, with unique ground states on a torus and 
no topological order. 

In this paper, we want to address the following ques- 
tions: Is the string order stable against small perturba- 
tions? And can it be generalized from one- to higher- 
dimensional systems? More precisely, we will investi- 
gate weakly coupled two- or three-dimensional arrays of 
one-dimensional systems. In the absence of the higher- 
dimensional coupling A_L they are characterized by vari- 
ous types of string order and a finite gap in the spectrum 
(but no topological order). Due to this gap, a sufficiently 
small Aj_ will never induce any thermodynamic phase 
transition. Therefore, one might naively expect that also 
the order parameter is robust against such a small pertur- 
bation. While this is correct for local order parameters, 
we will show that generically it does not hold for the 
non-local string order. 

In the following, we will first consider, as a prototypical 
example, the disordered phase of the quantum Ising chain 
in transverse field characterized by a hidden string order 
parameter. Here the divergence of the perturbation the- 
ory indicates that string order is destroyed by arbitrary 
small higher-dimensional coupling A_l. We connect this 
result - that also holds for the case of the spin-1 Haldane 
chain - to the band-insulators case where an exact calcu- 
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lation of the string order-parameter is possible, proving 
the absence of string order for any finite A_l. Using the 
language of bosonization, we identify the general mech- 
anism destabilizing non-local order. In contrast, local 
order (e.g. a charge density wave) remains stable. 

Transverse field Ising model: As a starting point for 
our discussion, we introduce the Hamiltonian of the 
quantum Ising chain 



H 



(3) 



describing a quantum magnet in a transverse field. We 
will always consider T = and the ferromagnetic case 
J > ( J < leads nevertheless to completely equiva- 
lent physics). The transverse-field Ising chain is a text- 
book example in the context of quantum criticality. Upon 
increasing the tuning parameter B/J, the ground state 
experiences a quantum phase transition from a magnetic 
to a paramagnetic state and the model can be solved ex- 
actly with the use of the standard Fcrmionic representa- 
tion for a quantum spiniS. Interestingly, there is a subtle 
way to identify the critical value oi B / J that exploits a 
hidden non-local property. In perfect analogy with the 
Kramers- Wannier^^ duality transformation for the clas- 
sical two-dimensional Ising model, one can introduce the 
following mappingi^ 



Mr 



(4) 



that preserves the SU (2) algebra and transforms Eq. ([3]) 
into 



H = Y,B^,l^,l 



J/i^ 



(5) 



i.e. one obtains the same Hamiltonian as in ([3]) when re- 
placing J <r^ B, ^ af . Therefore the quantum critical 
point has to be located at the self-dual point J/B = 1. 

As in the ferromagnetic phase B < J, (cr^) is finite, 
one finds in the disordered phase, B > J, that (/i^) is 
finite, or more precisely 



lim = lim 




>0 (6) 



Therefore the disordered phase is characterized by non- 
local string order 

If we now take a pair of these quantum Ising chains, 
we can investigate the fate of the non-local order when a 
weak inter-chain coupling J± is present. For simplicity, 
we set J — and consider the Hamiltonian 



H = J2Bal,+Bal, + J^al 



'2A- 



(7) 



As we have set J = 0, the system is a sum of independent 
two-site Hamiltonians and is trivial to calculate the string 



order exactly 



n ' 

yi<m<_; 



rli-«l 
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where we used that (cr^ 



J\/{%B'^) to leading order in J^_/B. The string order, 
being a product of factors strictly less then 1, decays ex- 
ponentially for any Jl 0. 

So far, the vanishing of the string order can be an ar- 
tifact due to the absence of interactions. To see that a 
finite J < B cannot stabilized the order, it is instruc- 
tive to repeat the calculation in the dual variables. In 
this language, the order parameter is now local but the 
coupling Jj^ induces a non-local term in the Hamiltonian 



H = 



rn<i 



k<i 



(9) 

With standard text-book techniques, we expand the S- 
matrix to second order in J± and obtain for the two-point 
correlation function 



.71 



8B^ 



(10) 



consistent with ([8]). Formally, the divergence with |j — j| 
arises from the term 



Jo J-oo 7. , 



k.l m<ik 



S<1 

(11) 

as a consequence of the non-local commutation relation 
between the order parameter and the perturbation 



7n<l 



m<.l,m^i 

(12) 

Physically, nm<; A^m describes the creation of a domain 
wall. As any domain wall created between the points i 
and j destroys the correlations of yitf ^ and /xf ^ , the per- 
turbation theory diverges linearly in |i— j|. From the per- 
turbative argument, it is easy to see that a finite J < B 
does not change this picture qualitatively for \i — j\ large 
compared to the correlation length ^ ~ (S — J)~^. In 
fact, the interaction term ^ J(J^cTi_^i = '^A'f has local 
commutation relations with the order parameter and pro- 
duces only regular corrections to the perturbation series. 
These additional non-singular terms cannot compensate 
for the linear divergence induced by the inter-chain cou- 
pling. 

Even though the long range order vanishes, thermody- 
namics is unaffected by Jj_ (up to a small renormaliza- 
tion of the gap): obviously no quantum phase transition 
is induced. Only the order parameter is sensitive to the 
non-locality of the perturbations. In the language of the 
original variables af, in contrast, the perturbation was 
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local, the order parameter non-local, and the same type 
of divergence arises again due to their non-local commu- 
tation relations ([T2|) . 

Finally, notice that in the ferromagnetic phase, J > B, 
the local magnetic order is stabilized rather than sup- 
pressed by J±. While for B > J, J± induces a finite 
density of virtual 'dual' domain wall fluctuations, the 
'physical' domain walls are suppressed for B < J hy J±. 
More precisely, they are confined as the energy of a pair 
of domain walls with separation |j — j| is proportional to 
\i — A in the ferromagnetic phasei^. 

Spin 1 chain: A second major type of non-local or- 
der that we want to consider is the so-called string or- 
der. This was introduced first in 1989 by Den Nijs and 
RommelseS briefly after that Haldane conjectured the ex- 
istence of a gap in the spin-1 antiferromagnetic chain as 
the hidden order of the Haldane phase. They observed 
that, even though true Neel order is absent, the ground 
state has still a form of long range order characterizing 
the entanglement of the spins: any site with S*^ = ±1 
is always followed by another with = ^1, separated 
from the first by a string of S'^ = of arbitrary length. 
This implies that the string order parameter 

50chai„(* - j) - [S^ exp(z7r ^ 5f ) S^) (13) 

l=i+l 

is always finite for |z — j| — > cx) in the Haldane phase. 

Kennedy and Tasak i^^i^^ observed that this non-local 
order can be understood from the non-local mapping 

L 

k=j+i 

Sy = exp(z^^5,^)5jexp(z^ ^ S^,) 

k=l k=j + l 

^1 = exp(z^^^,^)5; (14) 

as it maps the Heisenberg Hamiltonian with open bound- 
ary conditions to an effective local ferromagnetic Hamil- 
tonian with a manifest Z2 x Z2 symmetry that is fully bro- 
ken in the ground-state (the origin of the 4-fold ground 
state degeneracy for such boundaries is well understood 
in the AKLT picturei^ and comes as a consequence of an 
effective spin 1/2 localized at each boundary). The same 
transformation, applied to the inter-chain coupling, in- 
troduces in the Hamiltonian domain-wall creation opera- 
tors similar to the ones of Eq. ^ . Therefore the absence 
of the string order in the presence of a finite Jl can be 
shown as above. Indeed, the vanishing of the string order 
for two coupled spin-1 chains was observed numerically 
by Todo ct alM- using the quantum Monte Carlo method. 

Ladders: Similar to the spin-1 chain, also gapped spin- 
1/2 ladders are characterized by a non-local string order—. 
In this case, one can distinguish between two different 



types of string orders, ^Oodd and S'Ooven (see Ref. [1]) 
^Oodd(*-.7) = -{iSl, + Sl,)exp{zn ^ S^^ 

l=i+l 

+Sh){Sl, + Sl^)) (15) 
50cvcn(i-j) = -((5r,,+i +5|^,)exp(i7r ^ Sli+^ 

l=i+l 

+Sli){Sl,+, + Sl^)). (16) 

where the spin-1 of Eq. (|13p is replaced by the sum of 
two spin-1/2 operators on either the vertical or on the 
diagonal rung. 

In a recent paper—, we have shown that the ground 
state of gapped spin-1/2 ladders (and of the spin-1 chain) 
is adiabatically connected to the one of an ordinary, non- 
interacting band insulator. In Ref. [3], we introduced the 
following family of ladder Hamiltonians 

+ X! ^«'^i,^cr«2,i,o- + i_DaI,j+i,o-a2,i,<T + h.C. 

+ C/ ^ na,i,]na,i^l -f Jr ^ Si,iS2,i (17) 

i,a i 

(where a — 1, 2 and a ~ ± are the row and spin indices) 
whose phase diagram includes both Mott {U,Jji ^ t^) 
and band insulating phases {U — Jr — 0). Remarkably, 
the string order (|15I16[) turns out to be finite also for the 
band insulator. 

In the non- interacting case {U = Jr = 0), the string 
order can be calculated exactly^ even for an array of such 
one-dimensional insulators. As shown in Fig. [Uin the 
presence of an arbitrarily weak inter-ladder coupling 

H±=tj_ ^ 0,2.i,l,a'^i-i.l + i-A + h.C (18) 
i.l,a,(7 

(where the extra index I labels the ladders), the string 
order decays exponentially 

S-O-e-^l'-J'l, aoctl (19) 

as in Eq. ([5]). The prefactor a is quadratic in the inter- 
ladder hopping t± (see inset of Fig. [1]) . 

Bosonization: To identify the precise origin of the 
fragility of string order, we now analyze a generic 
one-dimensional gapped phase using the language of 
(Abelian) bosonization. Within this approach, the rel- 
evant degrees of freedom are described by the fields <& 
and O obeying the non-local commutation relation 

[^{x),<d{x')\ = ie{x-x') (20) 
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FIG. 1: (Color online) The odd string order for a two- 
dimensional array of weakly coupled non-interacting band in- 
sulators decays as 50odd ~ e""'*"-". For t± = —0.5 we also 
show the exponential fit. Inset: q as a function of t\. Both 
plots are done for the set of parameters ti — 0, t2 = —0.6, 
tR = -1.5, to = -2 of the model defined in Eqns. (fTTlfTS)) . 



where 9{x) is the usual 0-function. Gapped phases typ- 
ically arise when one of the Bosonic fields, e.g. $, is 
locked by a 'relevant' non-linear interaction 



gcos(x$) (21) 



where K is the Luttinger liquid parameter and g and x 
parameterize the most relevant perturbation of the Lut- 
tinger liquid fixed point. The Hamiltonian has two sym- 
metries: First, it is invariant under the shift Q O + c 
by an arbitrary constant c reflecting the conservation of 
the 'charge' / dx^ (typically either the total or the to- 
tal number of electrons) . Second, the field $ is invariant 
under a shift -I- by a fixed number c^, reflect- 

ing charge quantization and giving rise to the constraint 
X = n27r/c0 with integer n. For example, in the case of 
a non-interacting band-insulator = x — V^, 

K = 1 and g is proportional to the 2fci? (where kp is the 
Fermi momentum) component of the periodic potential. 

In this language the origin of the string order is easy to 
understand. As the field $ is locked in one of the minima 
of the cosine term, any correlation function of the form 



O-rix - y) 



_ /gi7'I'(a:)g-i7*(y) 



(22) 



does not decay for |a: — yj 



and has a finite value 



lim 0^{x-y)^nlo. 

\x-y\^oo 



(23) 



For example, a 'string operator' as in Eq. (|15p is 
bosonized in the continuum limit by^ 



exp 



(^7 E SI 

l=i+l 



-s. 



2.1 



■ exp 



($,(a;) - $,(y)) 



(24) 



FIG. 2: (Color online) The string order parame- 
ter without including the boundary spins defined as 

= limj^oo ^exp (j.^ I^Ci + ^2,1) ) ^o'' ^ 
dimensional band insulator on a ladder (for the value of 
the parameters, see Fig. [1} as a function of 7. By con- 
struction fi^ — 1 for 7 = Ann. Inset: String order SO^ = 

hm,^^ {{Sl^ + Sl^) exp (i7ECi + ^l.i) (Sf.o + S|,o)) 
including the boundary spins. Notice that at 7 = 27r the 
S02tt is exactly zero while the has a small but finite 
value. This arises because the operator e^'^'^^' is independent 
of the spin configuration. 



where dx^s is proportional to the total spin density. Note 
that in Eqs. (fTB|) and (fT5|) besides the string opera- 
tor, also extra boundary spins at site i and j have been 
included in the definition of the string order. However, 
these terms are not essential, as discussed in the caption 
of Fig. [21 and string order is present as long as 0^{x — y) 
is finite for |x — y| — > 00. 

We can make use of this formalism to understand the 
consequence of an inter-chain coupling. In the presence of 
such a perturbation, the total 'charge' J dx^ on a single 
chain is not conserved any more. This is reflected in the 
appearance of the dual field Q in the Hamiltonian (more 
specifically, of the exponential e*^®^^)). For instance, for 
chains coupled by single-electron hopping one obtains 



-Ye^V^i 



*c.i±*.,i+ee,!±e,.,) 



X e 



27ra 

-iw^('i'c.i+i±*s,i+i+ec.i+i±es,i+i) 



(25) 



where the summation index I spans the different lad- 
ders and we introduced the spin and charge Bosonic 
fields (and their duals). In the last equality only the 
I L,iJr\,a component is shown. 
According to the commutation relations ((20|) , the oper- 
ator e*'^®(^) increases the $(2;') field by /3 for x' < x. In a 
semi-classical picture, e*'^®*^^^ therefore creates a domain 
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wall at X by shifting $ for x' < x from one minimum of 
the cosine to another (note that (3 is always an integer 
multiple of 27r/x as a consequence of charge quantiza- 
tion). 

In analogy to the calculation for the transverse field 
Ising model, Eq. pT|) . we can now proceed by calculat- 
ing the corrections to the string order parameter per- 
turbatively in Hj^ — J dxh±{x). Here we assume that 
the string order is calculated on the chain with index 
I — and is defined in terms of the spin field ^s,o, 
0^{x - y) ^ (e»7(*..o(^,o)-*e,o(?/,o))^_ j,^ second order 
we obtain 

1 (■+00 /•+00 

- / dxidt, / dx2dt2(re*''^*=''''"'°^"*='''^^'°"x 



Here, 



h±{xi,h)h^{x2,t2))^ (26) 
denotes the connected part of the correlation 



function (i.e. (e 



i7(*»,o(a:)-$s,o(y)) 



){hj_h±) has been sub- 



tracted) and the time-ordering T splits the integral in the 
four different contributions (^i < < t2), (^2 < < ti) 
{ti,t2 > 0) and {ti,t2 < 0). As our system is massive, 
correlations decay on scale 1/A, where A is the gap, and 
therefore most contributions to are only of order 
{t±/A)^. There is, however, one important exception: if 
a domain wall is created at time ti < and position xi 
with x 'ti xi ^ y and destroyed at time ^2 > 0, the order 
parameter changes by e^^'^^ as at time t — a. domain 
wall is enclosed between x and y. Here we assume that 
|a; — yl, \xi — x\ and \xi — y\ are much larger than the 
correlation length ^ ~ We therefore obtain (up to 

corrections of order ^/\x — y\) 



exponentially small. As the non-local string correlations 
are destroyed by domain walls, the string order vanishes 
exponentially, 



SO 



,,2 p-a\x- 



-v\ 



(28) 



where a can be read off from Eq. ()27ll in the limit of small 
t± < A. 

In Fig. [3l we compare Eq. (p7)) with the exact calcu- 
lation available for the weakly coupled band insulators 
{rj — 1, P — ^7r/2). We find nice agreement with the 
expected a w 1/A^ and a « 1/A behavior both deep 
in the gapped phase, where ^ « a (left upper panel of 
Fig. [3]) and close to the critical point where £, ^ a (right 
upper panel of Fig. [3]) . Also the angular dependence is 
very well described by the factor (1 — cos(7/3)) (see lower 
panel of Fig. [3]). 

For the case of the quantum Ising chain in transverse 
field, the J = case has already been considered in 
Eq. dS]). We can now use our result to predict the leading 
behavior of the decay exponent for i? > J. As the scal- 
ing dimension for the SzSz coupling between neighboring 
chains is giveni^ by 77 = 1/4, it follows from Eq. ((27l) that 
the string order decays as 



SO{x - y) 



(29) 



where c is a dimensionless constant of order 1. 

From Eq. ((27|) we can directly infer the conditions of 
stability of the string order, in fact, if the relation 



7/3 = 2?™ 



(30) 



1 

- J dxi dti dx2 di2(e±*''^*(^~^i)e=F*'"^*(s'-^2) - l)x 

^g.-y(*.,o(.)-*.,ote))^^ {h^(xuh)h^ix2,t2)) 

^ M5ol2^-y|(cos(7/3) - l)x 

dti / dt2 / dx' {h^{x',ti)h^{Q,t2)) 



00 ^0 J —00 

.,|„fe_J^.(o»h«-l) (1)^(1)'"" (27, 

where is the order parameter in absence of the per- 
turbation, a is the lattice spacing, z the number of neigh- 
boring chains and rj is the scaling dimension of h±. The 
last factor {£,/a)^~'^^ is only of relevance in the asymp- 
totic regime |a: — ?/| » f ^ a, i.e. for a gap A small 
compared to the band width. 

In the cases of the transverse-field Ising model with 
J — and the non-interacting band insulator, we have 
shown that the linear correction (P7|) resums to an expo- 
nential. Physically, it is clear that this also will happen 
for the arbitrary systems considered above: The finite tj_ 
induces a finite density of domain walls and the probabil- 
ity of having no domain walls between x and y is therefore 



holds, the order parameter commutes with h± and re- 
mains finite. For example, in the case of the band in- 
sulator, we have /3 = 7r/2 and the string order is 
formally stable for 7 — Ann in Eq. However, 

in this case the operator exp ^17 J^iZl+i ^ii + ^21) 
comes trivially the identity and the correlation function 
exp (i7ED+i + Sl^) 5fi) = (5|i5fi) decays 
exponentially even for the purely one-dimensional model 
(see Fig. [2]). Therefore no stable order parameter exists 
for a band insulator (or the Haldane chain). 

As a consistency check, we now analyze a case where 
the order is purely local and therefore stable with respect 
to small perturbations. For example, spinless Fermions 
form a charge density wave for sufficiently strong interac- 
tions (of sufficient long range) . Within bosonization such 
a system is described hy^ 



H 



dx 



K 



(31) 

and for if > 1/2 the cosine term becomes relevant lock- 
ing the $ field, which implies a spontaneous breaking 
of the translational invariance for the underlying lattice 
model. Equivalently, (|3T|) describes the physics of the 
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FIG. 3: (Color online) We test the different limits of Eq. pTjl 
in the case of an array of band insulators (the Hamiltonian 
is defined in Eqns. (|17ll8p ) where we can calculate the SO 
parameters exactly (50odd/even ~ e~"'*~^'). Upper left panel: 
a multiplied by the ratio /S^ /t\ versus 1/A in the limit of 
A much bigger then the band width ~ a) for the set of 
parameters ti — 0, t2 = 0, tn = —1, tn = [—80,-6], t± = 
—0.2. The corrections to Eq. (|27p are of order 1/A. Upper 
right panel: a multiplied by the ratio A/t\^ versus 1/A in 
the scaling limit (^ 3> a) for the set of parameters ti = 0, 
t2 ^ 0, tR ^ -1, to = [-1.2,-1.05], t± = -0.02. The 
corrections to Eq. (I27|l are of order A. Lower panel: a divided 
by the predicted angular dependence (1 — cos (7/?)) for the set 
of para mete rs ti = 0, t2 = 0, tR = —1, to = —4, t± — —0.4 
(/3 = \/tvJ2 is fixed by the low energy expression of the inter- 
chain coupling and we vary 7 in Eq. (I22|l '). Lower panel inset: 
the decay exponent versus 7/3 for the same set of parameters. 

XXZ spin-1/2 chain where the condition K < ^ trans- 
lates to Jz > Jxy In this case, long-ranged Neel order 
develops in the ground state. In both cases the rele- 
vant order parameter, the staggered component of the 
Fermionic density ^'^5'^, is local. Within bosonization 
it can be extracted from Um|^_j,|^oo(e'2v^(*(^)-*(y))). 
A perturbation due to hopping to the neighboring chain 
or due to a spin-flip is proportional to e^*"^^. This im- 
plies 7 = 2^/^^, (3 = y/w and therefore 7/? = 2tt, fulfilling 
the condition of stability (see Eq. ([30]) ): the charge den- 
sity wave or the Neel order are "true" local orders in 
a gapped system which are stable with respect to small 
perturbations. 

In conclusion, we have shown that the string order of 
gapped one-dimensional systems is very fragile: an arbi- 



trarily small coupling to neighboring chains or ladders is 
sufficient for its vanishing. This has to be contrasted with 
the behavior of essentially all other ground state proper- 
ties which are minimally affected. The finite energy gap 
protects them such that a critical coupling is needed to 
induce a quantum phase transition. 

One may think about possible generalizations of the 
string order parameter for higher dimensional systems. 
Indeed, for two coupled spin-1 chains, Todoii suggested 
a generalization of the string-order parameter (also used 
in Ref. [H] to characterize the phases of frustrated spin-1 
chains) 

S02^ hm 5i%52%e^'=-+i(^M+^'.2)5i%-52%- (32) 

where Sx.y of Eq. ^ now includes the spins of both 
chains. The stability of SO2 comes from the fact that 
all non-local commutation relations of the string-order 
and the inter-chain coupling vanish as ^ S^^ commutes 
with H± . In a completely analogous way, one can define 
for any finite number of chains N a generalized order pa- 
rameter SOn (with or without boundary spins) that is 
non-zero. However, a generalization of the above formula 
for a two- or three-dimensional system is not possible. 
First, SOn vanishes exponentially even in the absence of 
any coupling between the chains. Second, if the string of 
Eq. (12) is generalized to a square, a cube or any other 
finite subset of spins, this will immediately lead again 
to non-local commutation due the presence of 'dangling 
singlets' at the boundaries of such structures. As for di- 
mensions d > 1 the surface of any non-local structure of 
infinite extension is infinite, we conclude that no direct 
extensions of string order to higher dimensional systems 
can be deviced and generalized string orders will decay 
as e~"^ where A is the area of the boundary^*'. 

An interesting question for the future is the investi- 
gation of the stability of various types of topological or- 
der. For example, Senthil and Fisher^ have shown that 
the ground state degeneracy of the deconfined phase of 
a two-dimensional Z2 gauge theory is stable with re- 
spect to a small inter-layer coupling. Finally, a related 
problem is the stability of various types of entanglement 
measure a 1 . 
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